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Abstract: We discuss the interplay of neutrino oscillation and decay properties at neu- 
trino telescopes. Motivated by recent unparticle scenarios, which open the possibility of 
new neutrino decay modes over astrophysical distances, we perform a complete classifica- 
tion of possible decay schemes, and we illustrate how different scenarios can be identified. 
Moreover, we show that the sensitivity of neutrino telescopes to standard neutrino prop- 
erties, such as the mass hierarchy or Sqp, is greatly enhanced in specific decay scenarios. 
In particular we discuss the impact of an astrophysical neutrino detection on terrestrial 
experiments, such as on the mass hierarchy measurement at NO^A. For example, we find 
that the scenario where only v\ is stable can be uniquely identified against all the other 
decay schemes, and that in this case CP violation can be established (for large #13) by the 
combination of Double Chooz with the track-to-shower ratio at a neutrino telescope, even 
if the flavor composition at the source is unknown. Our statements are based on a complete 
analysis of all the present solar, atmospheric, reactor and accelerator neutrino data, as well 
as on realistic simulation of future terrestrial neutrino oscillation experiments. 
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1. Introduction 

Neutrino telescopes [1, 2, 3, 4] are sensitive to neutrinos with an average energy and traveled 
distance many orders of magnitude larger than present neutrino experiments, and provide 
therefore a completely new window on both standard and non-standard neutrino properties. 
Apart from "conventional" parameters, such as neutrino masses and mixing angles, a very 
prominent example of such properties is the neutrino lifetime. Phenomenologically, from 
the observation of neutrinos from supernova 1987A, we know that at least one neutrino mass 
eigenstate must be stable over galactic distances. More stringent and explicit bounds can 
be derived from different observations when specific decay models are assumed (see, e.g., 
Refs. [5, 6, 7] for an overview). For example, solar neutrinos strongly limit the possibility 
of radiative decays [8], while for Majoron decays [9, 10] explicit bounds can be obtained 
from neutrinoless double-beta decay and supernovae [11]. Purely phenomenological (i.e., 
model-independent) bounds are, however, much weaker, leaving enough parameter space 
for the decay of any mass eigenstate over extragalactic distances [12, 13, 14]. Neutrino 
decays are usually described by a factor exp[— t/(rj)] = exp[— (Lmi) / (Eti)], where r, is 
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the rest frame lifetime of the i/j mass eigenstate boosted by 7 = E/rrii in the laboratory 
frame. Since L and E are given by the experiment and rrij is unknown, one typically quotes 
Ti/rrii as the neutrino lifetime. The best direct limit is obtained for v\ from SN1987A to 
be about T\jm\ > 10 5 s/eV [15]. Bounds on vi lifetime are imposed by solar neutrino 
data, yielding T2/rri2 > s/eV for invisible decays [12, 13] and T2/TO2 > 10~ 3 s/eV 

for decay modes with secondary v e appearance [16, 17]. Finally, limits on v-$ follow from 
the analysis of atmospheric and long-baseline neutrino data, T3/7713 > 10~ 10 s/eV [14]. 
Neutrino telescopes can probe lifetimes as long as ~ 10 7 s/eV • L [kpc] / E [TeV], which are 
about a factor of 10 5 longer than the current best limit (for E ~ 1 TeV). In view of the 
rather weak direct neutrino lifetime limits and the recently proposed unparticle models, 
which may lead to new mechanisms of neutrino decay [18, 19, 20, 21], we do not assume 
any specific decay model, but study the most general case. Note that neutrino telescopes 
may also probe different kinds of new physics (see, e.g., Refs. [7, 22, 23]), which we do not 
discuss in this work. 

In addition to decay properties, the propagation from the neutrino source to the detec- 
tor is described by the neutrino mixing parameters through averaged neutrino oscillations, 
i.e., by neutrino flavor mixing. If the neutrino telescope has some flavor identification 
capability, this dependence can be used to extract information on the decay [24, 25, 21] 
and oscillation [26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36] parameters, in a way which 
might be synergistic to terrestrial measurements. For example, reactor neutrino exper- 
iments in combination with astrophysical observations might provide hints on Sqp well 
before superbeams, by measuring the CP-even part of the oscillation probabilities [30, 32]. 
In addition, the sensitivity to the oscillation parameters can be enhanced in some decay 
scenarios [37]. In this study, we discuss the identification of the various neutrino properties 
in scenarios with both neutrino flavor mixing and decay. Since uncertainties in the mixing 
parameters [38] and flavor composition at the source [39] may limit such measurements, 
we carefully include these aspects in our study. 

Astrophysical sources and flavor composition. The existence of astrophysical 
neutrinos is not yet proven, but the detection of very high energy cosmic rays points towards 
cosmic accelerators which are expected to produce in addition high energy neutrinos. There 
are many potential candidates for neutrino sources, such as gamma ray bursts, active 
galactic nuclei or starburst galaxies, the latter being unaffected by the Waxmann-Bahcall 
bound [40]. Astrophysical neutrinos are normally assumed to originate from pion decays, 
with a flavor ratio at the source of (/ e , / M , f T ) ~ (1/3, 2/3, 0) arising from the decays of 
both primary pions and secondary muons ("pion beam source"); here f a is the fraction of 
flavor v a (neutrinos and antineutrinos combined), so that / e + /^ + /r = 1- However, it was 
pointed out in Ref. [41] that such sources may become opaque to muons at higher energies, 
in which case the flavor ratio at the source changes to (f e , f^, f T ) ~ (0, 1, 0) ("muon 
damped source"). Therefore, one can expect a smooth transition from one type of source 
to the other as a function of the neutrino energy [42, 43]. Once a specific neutrino source 
is found and identified, for example from its energy spectrum or using information from 
its optical counterpart, it might be possible to select a specific flavor ratio at the source by 
applying suitable energy cuts to the data. However, note that, since the neutrino flux drops 
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as the energy increases, we can expect less events from muon damped sources than from 
pion beam sources. In this work we will mainly focus on pion beam and on muon damped 
sources, as well as on neutrinos produced by photo-dissociation of heavy nuclei with a flavor 
composition (/ e , f^, f T ) ~ (1, 0, 0) at the source [44, 45] ("neutron beam source"). In some 
cases we will consider arbitrary flavor compositions at the source without significant tau 
neutrino production, as is expected from a diffuse flux coming from the superposition of pion 
beam and muon damped sources with different energy dependencies. For the identification 
of even more generalized sources including the possibility of tau neutrinos, see e.g. Ref. [46]. 

Detector and observables. On the detector side, flavor identification is the prereq- 
uisite to learn about neutrino properties. In a neutrino telescope such as IceCube, muon 
tracks can be most easily seen for E > 100 GeV. Electron and tau neutrinos will produce 
showers with a somewhat higher energy threshold, E > 1 TeV. In general, it is not possi- 
ble to distinguish between electron and tau events close to the threshold, whereas at much 
higher energies one may be able to identify these flavors as well [47]. In particular, one may 
use the "double-bang" signature of v T in a window 5-10 14 eV < E < 2-10 16 eV to distinguish 
all flavors [48]. It is therefore plausible to assume that the ratios R = 4>^ et / '((j)^ et + (j)^ et ) 
(tracks/showers) and S = 0g 6t /^ et (electromagnetic/hadronic showers) can be used as 
observables, where </>^ et = + d< ? is the neutrino (+) plus antineutrino (-) flux of flavor 
v a at the detector [30, 32]. As an additional observable, one may use T = </>g?/0^ ct for the 
Glashow resonance process v e + e~ — > W~ — > anything at around 6.3 PeV [48, 49, 29] to 
distinguish between neutrinos and antineutrinos. Note that T is detectable only in a very 
narrow energy range, and that due to the difficulties in the identification of the double-bang 
the precision on S will be lower than the one on R. Moreover, for muon damped sources 
S will only be measurable in rare cases, because the typical energy window of this kind of 
source coincides only occasionally with the double-bang window. 

This work is organized as follows. In Sec. 2 we describe the considered physics scenarios 
and we present a general classification of all possible decay schemes. In Sec. 3 we focus on 
normal mass hierarchy and discuss the possibility of identifying the decay scenario. The 
impact of the neutrino mass hierarchy and the possibility to establish it from astrophysical 
pion damped sources is discussed in Sec. 4. In Sec. 5 we extend our results to the case 
of unknown flavor compositions at the source, as is the case for certain diffuse fluxes. In 
Sec. 6 we study the Glashow resonance process as an additional observable, addressing 
the problem of the separation between neutrino and antineutrinos. In Sec. 7 we illustrate 
synergies with terrestrial neutrino oscillation experiments. Finally, in Sec. 8 we summarize 
our results and draw conclusions. Details of our statistics treatment can be found in 
Appendix A. 

2. Considered physics scenarios 

In this work we consider the most general combination between neutrino flavor mixing and 
arbitrary neutrino decay scenarios. This includes the conventional picture of flavor mixing 
among stable states as a limiting case. Concerning the oscillation part, we assume that 
flavor mixing take place only among the three known neutrino flavors, which means that 
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Table 1: Classification of all the possible decay scenarios for complete decays, according to both 
LMH (rows) and 123 (columns) naming conventions. The tags "L", "M", "H" refer to the lightest, 
middle, heaviest active mass eigenstate, respectively, whereas "I" refers to an invisible state. The 
tags "1", "2", "3" refer to the v\, V2, 1^3 active mass eigenstates. A slash through the tag means 
that the corresponding state is unstable. The icons illustrate the correspondence between the two 
naming conventions according to the given mass hierarchy. The black and white disks correspond 
to stable and unstable mass eigenstates, respectively. 

sterile neutrino states - if they exist - do not mix with the active ones. As for the decay 
part, following the approach of Ref. [24] we assume that all unstable mass eigenstates have 
decayed between the source and the detector, i.e., the decays are complete. Moreover, 
we neglect possible differences between neutrino and antineutrino decay rates, and we 
assume that if one polarity has completely decayed, the other one has as well. The decays 
products may be visible to the neutrino detector, i.e., different active states, or invisible 
for the detector, such as sterile neutrinos, unparticle states, Majorons etc.. Since the 
decay is assumed to be complete and neutrino oscillations are completely averaged over 
astrophysical distances, the transition probabilities are independent of the neutrino energy. 
Therefore, in this study we do not take into account possible information on the energy 
spectrum, and only focus on total rates. In particular, we assume that the daughter 
neutrinos, if active states, fully contribute to the observed signal regardless of whether 
they are degraded in energy. Note that in some scenarios interference effects between 
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oscillations and decay may occur if the source is coherent and the neutrinos decay while 
they are still oscillating [50, 51]; however, for the sake of simplicity we do not consider such 
cases or the corresponding corrections. 

Let us now consider all possible decay scenarios in a systematic way. First of all, 
note that for kinematical reasons any mass eigenstate can only decay into lighter ones, 
which in turn may be stable or unstable. However, the assumption of complete decay 
allows to eliminate intermediate unstable states from every decay chain. For example, if 
the heaviest eigenstate decays into both the middle and the lightest state, and the middle 
state decays into the lightest state, finally everything will end up in the lightest state. 
This argumentation includes more complicated scenarios with arbitrary branching ratios, 
including active states decaying into invisible states which then decay back into active ones, 
as long as the initial states are active. It means that the transition probabilities can be 
written in terms of the effective branching ratios Br^j between the initial unstable active 
states Vi and the final stable active states vf. 

P°f> = E (l^/| 2 + E l^| 2 Br W j \Up f \ 2 . (2.1) 

/ stable \ i unstable / 

Note that Ylf ^ T i-*f = 1 onr y if there are no invisible final states, whereas in general 
J^Br^/ < 1. Thus this formula also accounts for invisible final states. 

As we have seen, any decay scenario is uniquely characterized by the stability of its 
active states. In general, there are 2 3 = 8 possibilities, since either active state may be 
stable or not. We list these possibilities in Table 1, together with the relevant parameters 
needed to completely describe each scenario. It is convenient to classify the various decay 
scenarios according to two different naming conventions: 

LMH classification. This naming convention is illustrated in the rows of Table 1. The 
labels "L", "M", "H" refer to the lightest (vl), middle (vm), heaviest (uh) active 
mass eigenstate, respectively; a slash through the label {e.g., "^1", "|4") means 
that the corresponding state is unstable. Different scenarios are denoted by #n, with 
n = 1...8. 

123 classification. This naming convention corresponds to the columns of Table 1. The 
labels "1", "2", "3" refer to the mass eigenstates v±, 1^2, ^3 relevant for neutrino flavor 
mixing, irrespective of their mass ordering; as before, a slash through the label (e.g., 
"y$") means that the corresponding state is unstable. Different scenarios 
are denoted by $n, with n = 1 . . . 8. 

The correspondence between these two classifications depends on the neutrino mass hier- 
archy; it is illustrated in Table 1. Specifically, for the normal hierarchy we simply have 
#n = $n, whereas for the inverted hierarchy we have #1 = $1, #2 = $3, #3 = $4, 
#4 = $2, #5 = $6, #6 = $7, #7 = $5 and #8 = $8. Note that the branching ratios Brj^ F 
are completely independent from the neutrino mass hierarchy when / and F are written 
in the LMH notation. Similarly, the matrix elements | ^q,^ | 2 and |c7 a j| 2 , which appear in 
Eq. (2.1), are insensitive to the mass hierarchy when i and j are 123 tags. On the other 
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hand, combinations of both mixing angles and branching ratios, such as the transition prob- 
abilities P a/ 3 in Eq. (2.1), do depend on the mass hierarchy. We will discuss the impact of 
the neutrino mass hierarchy in greater detail in Sec. 4. 

The observables R and S defined in the previous section can be computed as 

p _ fa Potji q _ X/q fa Pge ,n \ 

u ~ T f (P + P V T f P ' { ' 

/ jn, J a K 1 ae ~ 1 ar ) L^ia J a 1 ar 

Note that the overall flux normalization, which depends on the source luminosity, distance 
to the source, etc., cancels in this definition. In addition, although the cross sections are 
not very well known at high energies, the valence quark contribution becomes negligible, 
and the dependence on the flavor becomes small [52, 53]. Therefore, we expect that the 
uncorrelated cross section error among the different flavors are small, whereas the correlated 
cross section error cancels in Eq. (2.2). 

An interesting special case is when there is only a single active stable neutrino mass 
eigenstate Vf. Substituting Eq. (2.1) in Eq. (2.2), it is straightforward to see that (see also 
Ref. [54]) 

R -\U ef \i + \U Tf \i-l-\U» f \*> S -\U Tf \> (2 ' 3) 
since the probabilities factorize in source-dependent and detector-dependent parts. This 
means that in scenarios #6 and #7 the quantities R and S do not depend on the parameters 
a or b listed in Table 1, even if the probabilities do. In addition, Eq. (2.3) implies that there 
is no dependence on the f a characterizing the source, which means that the uncertainties 
on the source flavor composition are irrelevant for scenarios #5, #6 and #7. 

To explicitly illustrate the kind of implications that the observation of an astrophysical 
neutrino source could have for neutrino phenomenology, we show in Fig. 1 the dependence 
of R on the parameters Sqp (for large sin 2 26*13) an d sin 2 6*23 (for small sin 2 2^13). For 
definiteness we focus on the normal mass hierarchy and a pion beam source. As soon as 
a particular physics scenario is identified, a concrete measurement of R may considerably 
help in the determination of 5cp or 6*23. In turn, if 5cp, #13, and 6*23 are constrained, we 
can use R to infer the decay scenario. For what concerns 5qp, scenarios #1, #2 and #5 are 
obviously unfortunate, since there is little or no dependence on this parameter. However, 
if nature has implemented scenarios #3, #4, #6 or #7 astrophysical sources may provide 
very important information on 5cp- As for #23, and in particular the octant determination, 
all scenarios may help, but scenarios #4, #5, and #6 are especially well suited. Note that 
there is some parameter dependence on the branching ratios in scenarios #2 and #3, while 
there is no dependence on the flavor composition at the source for scenarios #5, #6, and 
#7 as explained above. Such synergies between astrophysical and terrestrial experiments 
will be discussed in detail in the next sections. 

Sometimes it may be useful to test additional assumptions coming from specific models. 
In addition to the general case, in this work we will consider the following special cases: 

Special case 1: the lightest mass eigenstate ul is stable. This constraint might be mo- 
tivated by the observation of neutrinos from supernova 1987A. Only scenarios #1, 
#2, #3, and #7 are compatible with this assumption. 
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Special case 2: there are no invisible states. It follows that the lightest state must be 
stable, since it could only decay into a sterile state. Therefore, this case is a special 
realization of the previous one. In addition, the branching parameters become con- 
strained: scenario #2 has now only one parameter (b = 1 — a), and scenario #3 has 
no parameters at all (a = 1). 

3. Physics scenario identification 

Let us now focus on the normal mass hierarchy (we will discuss the impact of the mass 
hierarchy in the next section) and let us assume that we can identify the source. Further- 
more, we assume that information on both observables R (muon tracks to showers) and 
S (electromagnetic to hadronic showers) will be available. Under these hypotheses, the 
measurement of an astrophysical neutrino flux corresponds to a point in the (R, S) plane, 
with certain measurement errors. 

In order to discuss to which extent one can in principle disentangle different physics 
scenarios, we show in Fig. 2 the 99% allowed regions corresponding to different decay 
scenarios. Specifically, we project the global x 2 from present and future terrestrial experi- 
ments onto the (R, S) plane for each scenario, as the oscillation parameters and branching 
ratio parameters a and b (where applicable) are varied. The left and right columns cor- 
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R R 

Figure 2: Allowed regions at 99% CL in the (R, S) plane corresponding to different decay scenarios, 
for a muon damped source (left panels) and a pion beam source (right panels). We assume a normal 
hierarchy. The upper panels correspond to the analysis of present data reported in Ref . [55] . The 
other panels show the impact of 3 years of Double Chooz data taking (1.5 with near detector), 
assuming no signal (sin 2 26 , 13 = 0, middle panels) or a large signal (sin 2 26*13 = 0.1, lower panels). 
The extra branching ratio parameters a and b have been varied as well, where applicable. 
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respond to muon damped and pion beam sources, respectively. In the upper panels we 
show 99% regions implied by the global analysis of present solar, atmospheric, reactor and 
accelerator neutrino data [55]. These experiments are further combined with the accurate 
measurement of sin 2 29\s expected after 3 years of Double Chooz data taking (1.5 of these 
with near detector), assuming that no signal (sin 2 29\3 = 0, middle panels) or a large signal 
(sin 2 2#i3 = 0.1, lower panels) is observed. 

Let us first of all focus on a pion beam source and current experiments, i.e., the upper 
right panel in Fig. 2. Ignoring the uncertainties in the astrophysical measurement, scenarios 
#4, #5 and #7 are clearly separated from each other and from the rest. Although there is 
some overlap among the other scenarios, the physics can still be clearly identified in many 
case. For example, a measurement (R, S) ~ (1, 1.5) would uniquely determine scenario 
#6. Conversely, if one measures (R, S) ~ (0.6, 1) there is some ambiguity among scenarios 
#1, #3, and ^6, but scenarios #2, #4, #5, and #7 can still be excluded. Even if the 
observable S won't be available, in many cases one can clearly identify or exclude certain 
scenarios using only the projection onto the i?-axis. For example, a measurement R = 1.5 
can only arise from scenario #5. The situation obviously improves if one includes in the 
analysis future terrestrial experiments such as Double Chooz, since the regions become 
somewhat smaller, although the qualitative picture does not change. Note, however, that 
the standard oscillation scenario #1 can never be uniquely established from astrophysical 
measurements, unless a further hypothesis on the stability of certain mass eigenstates are 
assumed a priori. 

If we compare the pion beam source (upper right panel) with a muon damped source 
(upper left panel), we see that in general the allowed regions are much larger for the 
muon damped source. This means that the physics scenario identification becomes more 
difficult, but it also implies that the dependence on the individual parameters is stronger. 
Indeed it is well known that in the standard oscillation case (#1) the dependence of R 
and S on sin 2 #23 and 5qp is considerably stronger for the muon damped source than 
for the pion beam source. It is also clear that the information from different sources is 
somewhat synergistic for what concerns the physics scenario identification. For example, 
(R, S) = (0.7, 0.5) single out scenario #4 for a pion beam source, whereas the scenario 
cannot be determined for a muon beam source. In turn, (i?, S) = (1.2, 0.2) not only points 
towards scenario #3, but also uniquely identifies the source as muon damped, since for 
a pion beam source no region is present in this point of the parameter space. Moreover, 
note that scenarios #5, #6 and #7, which only have one stable active mass eigenstate, are 
independent of the source type, as we have pointed out in Eq. (2.3). 

As the next step, let us simulate a realistic astrophysical measurement with concrete 
statistics. The details of our simple simulation are given in Appendix A. Note that we 
normalize the source luminosity to the number L of muon tracks which would be observed 
in the detector in the absence of neutrino decay (scenario #1) and for sin 2 26*13 = 0. This 
means that we are properly using the same source luminosity for all the different scenarios. 
In general, a la error of order 10% might be expected for (9(100) events [30]. Since this 
error is much smaller than the typical size of the allowed regions shown in Fig. 2, we 
expect that our considerations hold as long as there are enough events. This statement is 
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Table 2: Exclusion of different physics scenarios from present data plus astrophysical measure- 
ments, assuming normal mass hierarchy and a pion beam source. The first three columns refer to 
the simulated scenario and the simulated (benchmark) R and S values corresponding to the current 
best-fit values for the oscillation parameters and a = b = 0.5 (where applicable). These points are 
marked in the upper right panel of Fig. 2. The next seven columns give the Ay 2 at which the 
corresponding decay scenario can be excluded, marginalized over all branching ratios and oscilla- 
tion parameters. In the last two columns we marginalize also with respect to the fitted "wrong" 
scenarios. The different groups assume different statistics (muon tracks L) and either R and S as 
observables, or only R (c/., Appendix A). The % 2 from present solar, atmospheric, reactor, and 
accelerator data has been added [55]. 

quantified in Table 2, which is based on the combination of present data with a simulated 
astrophysical measurement. The first three columns refer to the simulated scenario and 
the simulated R and S values, which are also plotted in the upper right panel of Fig. 2. 
The next seven columns give the Ay; 2 at which the corresponding decay scenario can be 
excluded, marginalized over all branching ratios and oscillation parameters. In the last 
two columns we marginalize also with respect to the fitted scenarios, giving the Ax 2 and 
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number of sigmas at which the simulated scenario can be established against all the others. 
The various blocks refer to different assumptions about statistics (muon tracks L) and 
observables (either R alone or both R and S measured) . This table clearly illustrates that 
even for low statistics (middle block) many scenarios can be excluded. For example, for 
simulated scenario #1 and only 10 muon tracks, scenarios #3 and #7 are ruled out at 
more than 3a. The same conclusion holds for higher statistics but no S measurement 
(lower block). In the high statistics case with both R and S measured (upper block) it is 
possible to uniquely establish scenarios #4, #5, and #7, excluding all the other ones. 

An interesting and somewhat simpler issue is whether we can establish the stability 
of the lightest neutrino, i.e., if special case 1 of Sec. 2 is realized. This corresponds to 
the following questions: if the true scenario is one of #1, #2, #3, and #7 (in which the 
lightest neutrino is stable), can we rule out scenarios #4, $5, and #6? Conversely, if the 
real scenario is one of #4, #5, and #6 (in which the lightest neutrino is unstable), can 
we rule out scenarios #1, #2, #3, and #7? Focusing again on normal hierarchy and on a 
pion beam source, and assuming that both R and S can be measured, we can answer these 
questions by looking at the upper right panel of Fig. 2. Qualitatively, the perspectives to 
establish the stability of the lightest neutrino state depend on the true decay scenario, as 
follows: 

#1: never (it is contained in #6); 

#2: always (apart from a very small overlap with #6), if statistics is good enough; 
#3: sometimes, if no overlap with #6 and there is enough statistics; 
#4: always, if statistics is good enough; 
#5: always, even with low statistics; 

#6: sometimes, if no overlap with #2 and #3 and there is enough statistics; 
#7: always, even with low statistics. 

Therefore, the chances to determine whether the lightest state is stable or not over astro- 
physical distances are quite high. 

4. Neutrino mass hierarchy 

As we have seen in Sec. 2, the neutrino mass hierarchy plays a crucial role in the correspon- 
dence between the two different naming conventions introduced in Table 1. Therefore, a 
detailed discussion of the main features of each scheme can help to understand the impact 
of the mass hierarchy on the scenario identification discussed so far. In the language of the 
LMH scheme, the branching ratios - and therefore the propagation from the astrophysical 
source to the neutrino detector - are completely independent of the hierarchy. On the other 
hand, the projection of the flavor states (z/ e , v^, v T ) onto the mass eigenstates (vl, I'm, v h) 
is different between normal and inverted hierarchy. This means that the production and 
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Figure 3: Allowed regions at 99% CL in the (R, S) plane corresponding to different decay scenarios, 
for the normal hierarchy (left panels) and the inverted hierarchy (right panels). We assume a pion 
beam source. The upper panels show the general case where all the scenario are allowed, while the 
lower panels correspond to the assumption that lightest mass eigenstate is stable (Special Case 1, 
see Sec. 2). If we further strengthen this restriction by imposing that there are no invisible states 
(Special Case 2), the patterned regions can be excluded as well. 

detection processes, when described in terms of (vl, vm, ^h), look different in the two 
hierarchies. For example, for the normal hierarchy, u e is mostly vl, but for the inverted 
hierarchy, it is mostly vm- Therefore, the LMH convention provides a simpler description 
of the decay scenarios and is also more motivated from a theoretical (decay model) point 
of view, but it is not appropriate to describe oscillation phenomena. On the other hand, 
in the language of the 123 scheme, the relation between the flavor states (u e , ZAj, u t ) and 
the mass states (y\, U2, ^3) is the same for both hierarchies, but there is an asymmetry 
in the branching ratios created by kinematics, since heavier states can only decay into 
lighter ones. As an example, let us consider scenario $2 {v\ and vi stable, ^3 unstable). 
For the normal hierarchy, is the heaviest state and can decay into v\, into z^, or into 
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Table 3: Identification of the neutrino mass hierarchy from astrophysical measurements only. The 
first four columns refer to the simulated scenario, hierarchy and benchmark R and S values (marked 
in Fig. 3, lower row). The last six columns represent the overall Ay 2 and a for the wrong hierarchy 
exclusion marginalized over all physics scenarios. In these columns, we distinguish A%q and <7o 
for no special assumptions, A% 2 and o\ for special case 1 in Sec. 2 (lightest state stable), and 
A%2 and oi for special case 2 in Sec. 2 (no invisible states). We assume L = 100 muon tracks for 
this simulation of a pion beam source, and S and R to be measured. The \ 2 from present solar, 
atmospheric, reactor, and accelerator data has been added [55]. 

invisible states, with a plethora of branching possibilities which are described in terms of 
two parameters (cf. Table 1). Conversely, for the inverted hierarchy, 1/3 is the lightest state 
and can only decay into invisible states without any further freedom. This implies that 
for scenario $2 the branching possibilities for the inverted hierarchy are only a subset of 
those for the normal hierarchy, and hence the allowed region in the (R, S) plane for the 
inverted hierarchy is a subregion of the corresponding one for normal hierarchy. Similar 
analogies can be derived for scenarios $3 and $4. Note that scenario $1 does not depend on 
the branching ratios and therefore there is no asymmetry between the normal and inverted 
hierarchy regions. The same happens for scenarios $5, $6 and $7, since in this case the 
expressions for R and S are given by Eq. (2.3), which does not contain the branchings. 
Note, however, that there are very small discrepancies between the two hierarchies due to 
the slight asymmetry introduced by present data [55]. 

Let us now discuss the the impact of the neutrino mass hierarchy on the physics scenario 
identification. In Fig. 3 we show the allowed regions in the (R, S) plane corresponding to 
different decay scenarios, for the normal hierarchy (left panels) and the inverted hierarchy 
(right panels). The colors represent the different scenarios in the LMH classification, but 
each region is explicitly labeled according to both schemes. From this figure we observe that 
the allowed domains corresponding to the same 123 scenario but to different hierarchies 
are quite similar, and one of the two is always a subregion of the other. Therefore, all the 
considerations presented in Sec. 3 about the identification of the decay scenario in the case 
of normal mass hierarchy are still qualitatively valid for the case of unknown hierarchy, 
provided that they are reformulated in the language of the 123 scheme. In other words, 
by measuring the (R, S) parameters we may be able to uniquely establish which of the 
(u\, i>2, ^3) eigenstates are stable, but in order to convert this into a statement on the 
stability of (vl, um, vh) we need to know the mass hierarchy from an external source, 
such as long-baseline neutrino oscillation experiments (cf., Sec. 7.3). 

For what concerns the mass hierarchy determination, let us first consider the generic 
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case with all possible decay scenarios, shown in the upper panels of Fig. 3. The fact that in 
a given 123 scenario the allowed region for one hierarchy is always a subregion of the one 
for the other hierarchy adds to the already mentioned problem of the degeneracy between 
different scenarios, and hence there is only a very limited portion of the parameter space 
where the hierarchy can be determined unambiguously by an astrophysical measurement. 
On the other hand, at the end of Sec. 2 we discussed a number of special cases reducing the 
number of possible decay scenarios and also restricting the corresponding parameter space. 
These special cases can be either motivated by specific decay models, or by phenomeno- 
logical observations. For example, one may assume that the lightest mass eigenstate vl 
is stable (special case 1). This constraint implies that only scenarios #1, #2, #3, and 
7^7 are remaining, as shown in the lower row of Fig. 3. If one further assumes that there 
are no invisible states (special case 2), then the allowed branching ratios become more 
restricted and the patterned regions disappear as well. Clearly, the mass hierarchy can be 
now easily determined: for example, R > 1.2 would imply scenario #7 and the inverted 
hierarchy. Note, however, that this interpretation of the experimental result is no longer 
purely phenomenological, and is intrinsically linked to the special assumption used. 

We quantify this observation for several benchmark points in Table 3. The first four 
columns refer to the simulated scenario, hierarchy, and benchmark R and S values (marked 
in Fig. 3, lower panels). The last six columns represent the overall Ax 2 and a for the re- 
jection of the wrong hierarchy, marginalized over all physics scenarios. Here we distinguish 
A%q and do for no special assumptions, Ax? and o\ for special case 1 {vl stable), and 
A%2 and 02 for special case 2 (no invisible states). As can be seen, in the general case at 
most a 2a mass hierarchy determination is possible, even for the relatively high luminosity 
considered here. On the other hand, in special case 1 the mass hierarchy can be easily 
measured in most of the discussed cases. 

5. Generalized source or diffuse flux 

Now what happens if we do not know anything about the source, such as if we have a 
mixture of different sources, or even a diffuse flux? Can we still learn something about 
physics? Let us assume a flavor composition at source (X e , 1 — X e , 0), i.e., X e = f e is the 
electron (flavor) fraction, and there are no v^s produced. Such a flavor composition might 
be observed for a combination of different sources with different energy dependencies, or 
a diffuse flux. In these cases, X e can be obtained as a (weighted) average of the different 
X l e from the different sources i. In the most general case, we have < X e < 1, where 
X e = corresponds to a jU-damped source, X e ~ 1/3 to a pion beam source, and X e = 1 
to a source from neutron decays. Assuming that the neutrinos are only produced by pion 
decays (and partly subsequent muon decays) with an unknown energy dependence, we 
have < X e < 0.35 from Refs. [25, 39] including spectral effects. In general, any value of 
X e is possible, but only one physics scenario will be realized if the decays are complete. 
Furthermore, let us assume, for the sake of simplicity, that we know the mass hierarchy 
from a different source. 
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R R 

Figure 4: Allowed regions in the (R, S) plane for the normal hierarchy and an unknown source. 
The left panel corresponds to the electron fraction X e marginalized over in the range < X e < 1/3, 
while the right panel corresponds to it marginalized over the full range < X e < 1. See main text 
for details. 

We show in Fig. 4 the allowed regions for the observables R and S for the scenarios 
from Table 1 for the normal hierarchy (99% CL). Let us first of all assume that we do 
not know anything about the source(s), i.e., < X e < 1. Therefore, we marginalize in 
the right panel of Fig. 4 over X e in the full range < X e < 1, which means that the 
regions span the whole range between muon damped and neutron beam source. 1 As the 
first observation, the scenario with only one final active stable state remains unchanged in 
consistency with Eq. (2.3). As a consequence, scenario #5 is still easy to identify. For the 
rest of the scenarios there is relatively strong overlap, and only in rare cases the scenarios 
might be identified. Nevertheless, many scenarios can still be excluded. 

If we assume that only pion beam and muon damped sources (and mixtures of these) 
contribute, we find the result in the left panel of Fig. 4. Such an mixture might be mea- 
sured for very limited energy resolution, unknown source parameters of a specific source, 
or a diffuse flux in a certain energy range. In this case, the result is qualitatively not ex- 
tremely different from the previous discussion. For example, scenarios ^5 and #7 are still 
relatively easy to identify. In addition, the conclusions from the previous chapters remain 
qualitatively unchanged. For a quantitative update, see Table 4, which is similar to Table 2, 
but for different sources, and X e marginalized in the indicated ranges. That means that 
here the source is assumed to be unknown to some degree, or one measures a diffuse flux 

For arbitrary marginalizations < X e < X™ ax , and arbritrary fixed X e , see movies in Appendix B. 
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Table 4: Same as Table 2 for different sources, and X e marginalized in the indicated ranges. That 
means that here the source is assumed to be unknown to some degree, or one measures a diffuse flux 
(superposition of sources). Here L = 1 00 and a measurement of R and 5* is assumed. The different 
(simulated) benchmark points are marked in Fig. 4. The \ 2 from present solar, atmospheric, reactor, 
and accelerator data has been added [55]. 



(superposition of sources). Here L = 100, and a simultaneous measurement of R and S is 
assumed. The different (simulated) benchmark points are marked in Fig. 4, where the first 
group corresponds to the triangles A and squares □ the left panel (simulated X e = 1/3), 
the second group to the triangles V and squares □ in the left panel (simulated X e = 0), 
and the third group to the circles O in the right panel (simulated X e = l). 2 It is now 
very interesting to compare the first group in Table 4 to the first group in Table 2, which 
are only different by the marginalization over X e . While in some cases the result does not 
change at all (such as for simulated scenario #1 and fit scenario #5), the sensitivity is in 
some cases completely destroyed (such as for simulated scenario #4 and fit scenario #3). 
This can be easily understood from Fig. 4, since the corresponding regions now overlap 
each other. Similar results are obtained for the simulated muon damped source in the 
middle row of Table 4. For the neutron beam source in the last row of Table 4 there is, 

2 Boxes are used if both points concide, which is the case for all scenarios with only one stable mass 
eigenstate. In this case, the observables to not depend on the flavor composition at the source; c/., Eq. (2.3). 
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Figure 5: Allowed regions at 99% CL in the (T, S) plane, for a pion beam source and a normal 
hierarchy. The left panel corresponds to pp neutrino production, while the right panel corresponds 
to pry production. In the legend, the black and white disks correspond to stable and unstable mass 
eigenstates, respectively. 

however, a qualitative difference: Since the mass eigenstate v% is initially not populated 
for X e = 1 because we assume a simulated sin 2 26^ = (and therefore U e s = 0; c/., 
Eq. (2.1) and Eq. (2.2)), the stability of 1/3 is irrelevant, and models that differ only in that 
stability are physically equivalent. This means that the simulated models are paired, i.e., 
#1=#2, #3=#7, #4=#6, #5=#8 (and this last case is irrelevant since nothing arrives 
at the detector). 3 From Table 4, it is quite interesting that in scenario #3=#7, any other 
qualitative case can be significantly excluded even if the only assumption on the source is 
that there are almost no i/ T 's produced. 

6. Glashow resonance process as a third observable? 

The Glashow resonance process v e + e~ — > W~ — > anything at around 6.3 PeV [48] allows 
for the detection of electron antineutrinos only. Therefore, we define T = / a^ et as an 
additional observable. This is the only observable which is sensitive to the production of 
pions (and kaons) at the source by interactions of high energy protons with photons ("^7") 
or protons ("pp") [49]. In the jry process, mainly ir + are produced through the A reso- 
nance, which means that the flavor composition at the source is (/ e , /„, f T \ fg, f^, ff) ~ 

3 There can, however, be a difference in the fit \ 2 between two paired models, because we allow for 
sin 2 20i3 > in the fit. 
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(jp I 0, |, 0) (split up by neutrinos and antineutrinos with ^ /j = 1). In the process, 
a nearly equal mix between n + and 7r~ is produced, leading to (f e , f^, f T \ f s , f^, ff) ~ 
(|, |, 0||, ^, 0). If the detector is CP-blind and there is no asymmetry between neutrinos 
and antineutrinos, we can sum the flavor compositions of neutrinos and antineutrinos at the 
source in order to obtain (/ e , / M , f T ) = {f e + fg, U + fp,, fr + ff) - (3, §, 0) in both cases. 
Similarly, if the muons are damped at higher energies, we have (/ e , / M , f T \ f s , fp, ff) — 
(0, 1, | 0, 0, 0) for pj, and (/ e , /„, f T \ f g , f p , f f ) ~ (0, ±, | 0, ±, 0) for pp. In addition 
to the electron fraction X e describing the fraction of electron neutrinos (and electron an- 
tineutrinos), one can introduce a photon fraction X 7 describing the fraction of neutrinos 
produced by pj processes. In this case, the fraction 1 — Xj comes from pp interactions. If 
we assume that all electron neutrinos (antineutrinos) come from antimuon (muon) decays, 
one can parameterize the source as 



(/ej ffit It I fei f [ii fr) — 

1 + X~y 1 + Xy — X e — 3X^X e 

7, A e, o j U 



~-^X e , l-X-y-Xe + 3X,X e \ (6i) 

Conversely, the fraction of observed (useful) muon decays < £ < 1 is given by £ = 
X e /(l — 2X e ). Note that this parameterization can only describe the above decay chain, 
and is not useful for X e > 1/3, such as for a neutron beam source. 

We show in Fig. 5 the observables T and S for a pion beam source (X e = 1/3) and 
pp (left panel, X 7 = 0) versus pj interactions (right panel, X^ = l). 4 Note that neither R 
nor S depend on X 7 since the neutrino and antineutrino rates are added. Since the regions 
do not collapse to thin curves, i.e., T is a well-defined function T(S), there is obviously 
additional useful information in T. However, in the pp case, the different scenarios cluster 
along the diagonal and overlap each other in the same fashion as for S. Therefore, for 
scenario identification, T may not provide much new information. For the pj case, however, 
there is obviously new information. For example, if (T, S) = (0.1, 1) is measured, scenario 
#3 can be uniquely established, as well as the pj source can be identified. If only R and 
S were used, it would be fully contained in clusters #2 and #6 (c/., Fig. 2, upper right 
panel). Note that one can also establish the pp source in some cases. For example, if T > 4 
is observed, scenario #7 together with a pp source has to be realized. 



7. Synergies with terrestrial neutrino oscillation experiments 

In this section we discuss the synergies with the terrestrial neutrino oscillation experiments. 
We focus on two decay scenarios, which are, in wide ranges of the parameters, relatively 
easy to identify: scenario $5 (only 1/3 stable) and scenario $7 (only v\ stable). Since in 
both cases there is only one stable mass eigenstate, the observables depend on the mixing 
matrix elements only (c/., Eq. (2.3)). In particular, there is no dependence on the flavor 
composition at the source, branching ratios, and mass hierarchy. This means that there is 
also no mass hierarchy information from the astrophysical neutrino source(s). As discussed 

4 For arbitrary X 7 , see movies in Appendix B. 
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Figure 6: Comparison between physics scenarios $1 (no decay, left column), and $7 (y\ stable, right 
column), for Double Chooz plus an astrophysical source. Here a pion beam source L = 100 tracks 
as normalized luminosity for both sources is assumed, i.e., the left and right panels correspond 
to the same source luminosity. In addition, it is assumed that only R can be measured. The 
contours correspond to la, 2a, 3a (1 d.o.f.). The dashed curves are for fixing the other oscillation 
parameters. The current best-fit values and parameter errors are taken from Ref. [55]. The used 
simulated values are sin 2 2#i3 =0.1, Sqp — 7r/2, and a normal hierarchy. 

in Fig. 1, R varies strongly with Sep in scenario $7 (see also Ref. [37]), whereas R varies 
strongly with #23 in scenario $5. Therefore, we use these two scenarios in combination with 
terrestrial measurements, and compare them to the standard oscillation result scenario 
(no decays) or the same luminosity of the source. As for the time scale, we choose the 
next generations of reactor and long-baseline experiments. Namely we use Double Chooz, 
MINOS, and NOi/A as examples. In addition, we assume that there are no decay effects 
observed in terrestrial experiments, i.e., neutrino decay is visible only over astronomical 
distances. For details on the statistical simulation, see Appendix A. 

7.1 Can Double Chooz plus neutrino telescope measure <5cp? 

As in was pointed out in Refs. [30, 32], flavor ratio measurements might allow a measure- 
ment of 5qp already in combination with Double Chooz. However, for a pion beam source, 
which may be the most common one, the dependence of R and S on 5qp and the other 
oscillation parameters in the standard no-decay scenario is very moderate (c/., Fig. 2). In 
this case, knowledge from different sources, high statistics, and the use of different observ- 
ables is necessary to obtain useful information on 5qp- However, if neutrinos decay there 
can be a relatively strong dependence on Sep, depending on the specific scenario [37]. We 
demonstrate this effect quantitatively for a three year Double Chooz measurement and a 
relatively large sin 2 2#i3 in Fig. 6. In this figure, the precision in sin 2 2#i3-<5cp is shown 
for maximal CP violation implemented by nature. For the astrophysical source, we only 
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Figure 7: Left panel: Sensitivity to the #23 octant as a function of the true sin 2 2#i3 and true 
sin 2 023, where sensitivity is given at the 90% CL above/below the curves (as illustrated by the 
shadings). The different curves represent MINOS plus Double Chooz (DC) alone (thin gray curve), 
MINOS plus Double Chooz plus astrophysical information in scenario $1 (dashed curves), and 
MINOS plus Double Chooz plus astrophysical information in scenario $5 (thick solid curves). Right 
panel: Sensitivity to the normal mass hierarchy as a function of the true sin 2 2# 13 and true Sqp 
(stacked to the "Fraction of Sep"), where sensitvity is given at the 90% CL below the curves. The 
different curves represent NOz^A alone (thick gray curve), NOi/A plus astrophysical in scenario $1 
(thin black curve), NOz^A plus astrophysical in scenario $7 (thick black curve), and NOz/A plus 
astrophysical in scenario $7 with both R and S as observables (dashed curve). In both plots, 
L = 100 muon tracks were assumed for the flux normalization, only R was used as an observable 
(unless stated otherwise), and a normal hierarchy was simulated. The octant plot does not include 
the mixed (octant and sign) degeneracy. For details on the simulation, see Appendix A. 

assume a pion beam source producing 100 muon tracks for the standard scenario, and we 
only measure the observable R. The left and right panels correspond to the same source 
luminosity. Obviously, if neutrinos are stable, there will be hardly any information on 5cp- 
However, if only v\ is stable, even a 2a CP violation measurement might be possible (if 
the uncertainties on the other oscillation parameters can be further reduced, even 3a - see 
dashed curves). 

7.2 Octant determination with terrestrial experiments plus neutrino telescope 

We define sensitivity to the #23 octant if the wrong octant can be excluded at the 90% 
CL. Quite naturally, there will be no sensitivity if sin 2 #23 is close to 0.5, and there will be 
no sensitivity if the observables depend only on sin 2 2623 7 as it applies to the disappearane 
channel in long-baseline experiments. The improvement of the octant measurement using 
astrophysical neutrinos was discussed in Ref. [30] for no decays and a combination of 
terrestrial experiments. Here we focus on a shorter time scale. We assume that we have 
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information from Double Chooz and MINOS. We show in Fig. 7, left panel, the sensitivity 
to the 6*23 octant as a function of the true sin 2 26*13 and true sin 2 023, where sensitivity is 
given at the 90% CL above/below the curves (as illustrated by the shadings). The different 
curves represent MINOS plus Double Chooz alone (thin gray curve), MINOS plus Double 
Chooz plus astrophysical information in scenario $1 (dashed curves), and MINOS plus 
Double Chooz plus astrophysical information in scenario $5 (thick solid curves). We can 
read off this figure that there is no sensitivity to the octant in the 3a currently allowed 623 
range (marked in the figure). If R is measured for stable neutrinos ($1) and for a pion beam 
source, the potential substantially improves compared to no astrophysical information, and 
includes wide region of the currently allowed range. Typically, it depends on the true 
sin 2 2#i3. However, if only 1/3 is stable ($5), we can read off from Eq. (2.3) that 

sin 2 #23 cos 2 013 2 

R = 1 ~T~a TT~ ~ tan 023 ■ (7 - X) 

1 — sin #23 COS z V13 

This implies that 023 can be measured almost without parameter correlation by the neutrino 
telescope. In the left panel of Fig. 7 we show the excellent precision compared to scenario 
$1 for the same source luminosity and observable, where the sensitivity mainly comes from 
the astrophysical source. 



7.3 Mass hierarchy determination with NOi/A plus astrophysical 

For terrestrial long-baseline experiments, the mass hierarchy degeneracy [56], which deter- 
mines the mass hierarchy measurement, is, in general, located at a different value of c>cp 
than the original solution. In addition, it moves in the <5cp direction as a function of the 
true sin 2 20i3 (c/., Fig. 4 in Ref. [30] for NOvA). Since astrophysical neutrino sources are 
sensitive to cos<5cp, whereas first generation superbeams operated close to the oscillation 
maximum are mainly sensitive to sin Sep , the knowledge from the astrophysical source can 
improve the mass hierarchy measurement at the terrestrial experiments [30]. This effect is 
largest in scenarios where the dependence of the observables on <5cp is strongest. We illus- 
trate this behavior in the right panel of Fig. 7, where we show the sensitivity to the normal 
mass hierarchy as a function of the true sin 2 20i3 and true 5cp (stacked to the "Fraction 
of <5cp"). The different curves represent NOzaA alone (thick gray curve), NOi^A plus as- 
trophysical R in scenario $1 (thin black curve), NOi^A plus astrophysical R in scenario $7 
(thick black curve), and NOz/A plus astrophysical R and S in scenario $7 (dashed curve). 
We compare our results for scenarios $1 (no decays) and $7 (only 1^1 stable) for a pion beam 
source. Note that, as explained before, there is no sensitivity at all to the mass hierarchy 
from the astrophysical source alone. As can be seen from this figure, there is almost no 
improvement for stable neutrinos ($1), where the mass hierarchy can be determined for 
up to (depending in sin 2 20i3) 40% of all true 5qp- However, if only v\ is stable ($7), a 
measurement of R helps significantly, and the mass hierarchy can already be determined 
for up to about 65% of all true e>cp. The additional information on S, as shown as the 
dashed curve, helps even more. 
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8. Summary and conclusions 

In this study, we have discussed the identification of different decay and oscillation scenar- 
ios at neutrino telescopes. Furthermore, we have studied the measurement of the physics 
parameters within these scenarios by a neutrino telescope alone, and in combination with 
future terrestrial experiments. We have taken into account the present knowledge of the os- 
cillation parameters from a global fit of current solar, atmospheric, reactor, and accelerator 
data, we have statistically quantified the information from the astrophysical sources, and 
we have performed a complete simulation of future terrestrial experiments. For the observ- 
ables, we have mainly focused on the muon track to shower ratio and the electromagnetic 
to hadronic shower ratio, but we have also discussed the Glashow resonance process. We 
have performed a complete classification of effective decay scenarios and the corresponding 
branching ratios for complete decays, accounting also for possible invisible states. We have 
demonstrated that, depending on the physics scenario implemented by nature, the identi- 
fication of the scenario can be unique or ambiguous. For example, if only v\ or v% is stable 
(either of which can be the lightest depending on the hierarchy), the physics scenario can 
be easily identified. In the standard oscillation case, however, only specific scenarios can 
be excluded. 

As far as the impact of the mass hierarchy is concerned, we have demonstrated that 
one may be able to establish which of the (ui, zv 2 , 1/3) mass eigenstates are stable, but not 
their mass ordering, which determines the branching ratios (since only heavier states can 
decay into lighter ones). For example, without external mass hierarchy measurement such 
as from superbeams, one can in principle determine whether u\ is stable or not, but not 
if v\ is the lightest or middle mass eigenstate. This implies that a generic mass hierarchy 
identification is only possible in very small corners of the parameter space. However, 
if one imposes some model-dependent constraints then the mass hierarchy can be easily 
determined in most cases from astrophysical neutrinos alone. One possible such constraint 
is the assumption that the lightest neutrino mass eigenstate is stable. 

We have also studied the impact of flavor composition uncertainties at the source or 
the use of diffuse fluxes. For example, we have demonstrated that if no more than one 
active neutrino mass eigenstate is stable, there is no dependence of the observables on the 
flavor composition at the source. In order to study diffuse fluxes, we have marginalized 
the electron fraction at the source, which means that have taken into account arbitrary 
combinations of muon damped and pion beam sources. In this case, the physics scenario 
identification becomes quantitatively more difficult, but the qualitative conclusions still 
hold. Even if one allows for arbitrary production of v e and v„ neutrinos at the source, 
some physics scenarios can still be established. 

For what concerns the determination of the neutrino parameters in particular decay 
scenarios, we have chosen a number of examples in order to demonstrate the impact of an 
astrophysical measurement for future long-baseline and reactor experiments. For example, 
if a neutrino telescope measures the track to shower ratio from a pion beam source, then 
Double Chooz might be the first experiment to establish CP violation if v\ is the only stable 
state, whereas we have not found any CP violation sensitivity if all the neutrinos are stable. 
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As another example, we have demonstrated that there is some sensitivity to the #23 octant 
if MINOS and Double Chooz are combined with astrophysical data even if all neutrino mass 
eigenstates are stable, but if only 1/3 is stable there will be direct octant sensitivity from the 
astrophysical source alone. We have also illustrated how the mass hierarchy sensitivity at 
NOz^A would be enhanced by an astrophysical neutrino measurement if only v\ is stable. 
For large sin 2 2#i3, the fraction of <5cp, for which the hierarchy can be determined, could 
increase from 40% up to 80%. Note that this mass hierarchy determination is independent 
of any model-dependent assumptions on the decay scenarios. 

We conclude that an observation of astrophysical neutrinos at a neutrino telescope 
would be an important test of the oscillation and decay neutrino properties.. While it is 
difficult to obtain information on the neutrino lifetime without a distance measurement 
of the source, complete decay scenarios can in many cases be easily identified even if one 
takes into account the current measurement precisions of the oscillation parameters and 
uncertainties of the flavor composition at the source. Especially if neutrinos decay, the 
combination with terrestrial neutrino experiments may lead to early and surprising results 
even for the standard oscillation parameter measurements. An important prerequisite for 
such conclusions will be the flavor identification in the detector. 
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A. Statistical method and simulation 

A detailed description of our simulation of present solar, atmospheric, reactor and accel- 
erator neutrino experiment can be found in Ref. [55], from which we also take the cur- 
rent best-fit values and allowed parameter ranges. For the MINOS simulation, we follow 
Ref. [57] with a total luminosity of 5yr x 3.7 • 10 20 pot/yr and a 5.4 kt magnetized iron 
calorimeter [58] (the unit "pot/yr" refers to "protons on target per year"). For Double 
Chooz [59, 60], we use the simulation from Refs. [61, 62] with 1.5 years of data taking 
with far detector only, followed by 1.5 years with both detectors. For NOz^A, we use the 
simulation from Refs. [63, 30] updated to the numbers from Ref. [64] and a 15 kt detector 
mass. The future reactor and long-baseline experiments are simulated with the GLoBES 
software [65, 62]. 
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We define an astro-physical xLtro to be added to the GLoBES software or to present 
experiments as 



Xastro = nun 



2 E 



=1 



TiiO-Oi + Oiln 




(A.l) 



Here Tj corresponds to the theoretical (fit) rate and Oi to the observed (true) rate. The 
index i runs over all event types, such as muon tracks, showers, double-bang, etc.; £ is 
a source type-dependent unknown (free) flux normalization parameter to be marginalized 
over, and its error. Since we are not assuming any prior knowledge of the neutrino flux 
expected from the observed source, we conservatively set — > oo. 

We consider two different cases for T and O. If only R is measured we set n = 2 and 

Ti = ^ , T 2 =£ (iV e fit + iVf) , 

1 = A^ rue , 2 = iV e true + A^ rue , 

where Np is the (total) number of events for flavor up. If both R and S are measured we 
set n = 3 and 

Ti=£Af, T 2 = £iV e fit , T 3 = £iV T fit , 

Oi = iV^ ruc , 2 = iV e true , 3 = A^ ruc . 
The event rate for the flavor i/p in the detector is given as 

where / Q denotes the fraction of neutrinos produced as flavor v a at the source, e@ = eg/e^ 
is a relative efficiency compared to the muon track detection efficiency, 4> corresponds to a 
normalized luminosity at the detector, and (k) refers to the decay scenario in Table 1 (the 
probability is described by Eq. (2.1)). In order to compare different physics scenarios for 
the same source flux, we normalize to a number of muon tracks L observed in the detector 
for the standard oscillation scenario #1 and 6*13 = 0, i.e., 

(A.5) 



^ (1) 

a=l 



This normalization does not depend on the physics scenario. Therefore, it allows to compare 
different physics scenarios for the same source flux, and one can identify the physics in 
which one can most efficiently measure the target parameter. In addition, the number of 
observed muon tracks in the standard scenario is a quite intuitive one. For a flux close 
to the Waxmann-Bahcall bound, one may expect L = 0(100) muon tracks in about eight 
years [66]. 

For the sake of simplicity, we assume that e e = e T = 1. In a more realistic simulation 
one would probably have e e , e T <C 1, since the detector is sensitive to partially contained 
muon track events generated out of the fiducial volume, and the energy threshold for 
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\x events is lower [47]. Our assumption corresponds to choosing appropriate cuts such 
that e e ~ e T ~ 1, i.e., muon tracks and the other event types are detected with similar 
efficiencies. For the case of S, that of course implies relatively low event rates. In addition, 
we assume a background-free environment. Backgrounds could be easily included in our 
treatment, but they strongly depend on the source type, energy range, etc., whereas we 
want our simulation to be as much source-independent as possible. Finally, we neglect 
flavor identification uncertainties, whose impact on the determination of the parameter R 
has been estimated in Ref. [47] to be of the order of 20% - whereas for the parameters S 
and T it is considerably worse. 

Although our approximations are likely to have a non-negligible impact on our quan- 
titative results, it is quite easy to understand their implications at qualitative level. For 
example, the uncertainties in R and S due to flavor identification can be accounted for in 
Figs. 2-5 by describing the experimental results as extended ellipses rather than points: 
degeneracies among models will then arise when more than one region overlap with the 
experimental ellipse. Also, the main effects of backgrounds and systematic uncertainties is 
to reduce the significance of the fit, hence they are quite similar to a reduction in statistics, 
which we have explicitly illustrated in Table 2. Therefore, we believe that despite our 
simplications the main conclusions of our work are solid. 

B. Movies 

A number of movies corresponding to Sec. 2, Sec. 5, and Sec. 6 can be found at Ref. [67]. 

The first type of movies uses X e as a free parameter, i.e., the time parameter, and 
corresponds to Fig. 2. The flavor composition at the source is assumed to be (/ e , / M , f T ) = 
(X e , 1 — X e , 0). Note that X e = corresponds to a muon damped source, X e = 1/3 to 
a pion beam source, and X e = 1 to a neutron beam source. The movie shows the normal 
hierarchy (left) and the inverted hierarchy (right). Movie versions corresponding to all 
rows in Fig. 2 are available. 

The second type of movies is similar to the first type, but X e is marginalized over in 
the range < X e < A™ ax , and A™ ax is the time frame parameter. This movie corresponds 
to Fig. 4, but movie versions corresponding to the different rows in Fig. 2 are available. 

The third type of movies shows all observable pairs as a function of the photon fraction 
A 7 for a pion or muon damped source (versions for the different mass hierarchies are 
available as well). The photon fraction quantifies fraction of neutrinos is produced by j>7 
versus pp processes. For more details, see Sec. 6. 
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